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Abstract 

How impurity atoms move through a crystal is a fundamental and recurrent question in materials. The 
previous understanding of oxygen diffusion in titanium relied on interstitial lattice sites that were recently 
found to be unstable, making the diffusion pathways for oxygen unknown. Using first-principles quantum- 
mechanical methods, we find three oxygen interstitial sites in titanium, and quantify the multiple interpen- 
etrating networks for oxygen diffusion. Surprisingly, no single transition dominates, but all contribute to 
diffusion. 

PACS numbers: 61.72.S-, 66.30.J- 
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Controlling diffusion of impurity atoms into metals is a basic and ubiquitous technique of ma- 
terials design for millennia. The diffusion of oxygen in titanium impacts design of implant and 
aerospace alloys, as well as the formation of titanium-oxides. Increasing the oxygen content in 
titanium forms ordered layered-oxide phases, which rely on the diffusion of oxygen into alter- 
nating basal planes to form[[T]-[3l ; modeling the kinetics of ordering[4J needs information about 
diffusion. Initial stages of growth of titania nanotubes — e.g., for dye- sensitized solar cells — via 
anodization of a titanium metal substrate[5] involves the diffusion of oxygen. At moderate to high 
temperatures, oxygen diffuses into titanium from the surface oxide, making the metal brittle, and 
requires costly processing steps to avoid[6]. Designing titanium alloys with lower innate oxygen 
diffusivity has the potential to replace heavier alloys in aerospace to reduce greenhouse-gas emis- 
sions. Ultimately, understanding how to impede or accelerate the diffusion of oxygen requires a 
fundamental description of diffusion pathways through titanium. 

Diffusion of single oxygen atoms through hexagonal-closed packed (a) titanium initially ap- 
pears simple — proposed as atom-hopping between identical interstitial sites, following an Arrhe- 
nius relationship with temperature — but that simplicity hides a complex network of transition 
mechanisms. Oxygen prefers to occupy an octahedral interstitial site surrounded by six titanium 
atoms [7 J, and so modeling oxygen diffusion had assumed either direct octahedral-to-octahedral 
transitions through tetrahedral transition states (8]|, or from octahedral to metastable tetrahedral 
sites[9]. However, the recent discovery that the tetrahedral site is unstable in favor of a metastable 
hexahedral site lfTOl left an open question: how does interstitial oxygen diffuse through ar-titanium? 
Moreover, the ratio of oxygen diffusivity along basal (xy) directions and the c-axis (z) direction is 
nearly unity[8J despite no symmetry relationship between the basal plane and the c-axis. 

We explain these phenomena with first-principles predictions of three distinct interstitial oxy- 
gen sites in ar-titanium — including a new non-basal crowdion site with lower symmetry than the 
octahedral and hexahedral sites — that create three different lattices interconnected through a net- 
work of transitions. We derive the first entirely first-principles prediction of oxygen diffusion 
in a--titanium, and our predictions compare well with experiment. Our results explain the near 
isotropy of oxygen diffusion as transitions occur between the three different interstitial sites with 
similar absolute rates, so that nearly all portions of the transition network contribute to diffusion, 
rather than a single rate-controlling step. That is, all escape times contribute to diffusion, not just 
from the octahedral site. The metastable crowdion site and hexahedral site give a route for accel- 
erating oxygen diffusion. The full topology of the oxygen diffusion network highlights an unusual 
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complexity for the kinetics of interstitials in metals. 

The first-principles calculations are performed with VASP lfTTl[T2ll . a plane-wave density-functional 
theory (DFT) code. Ti and O are treated with ultrasoft Vanderbilt type pseudopotentials tfT^l [T41 
and the generalized gradient approximation of Perdew and Wang [fT5ll . We use a single oxygen 
atom in a 96-atom (4x4x3) titanium supercell with a 2 x 2 x 2 fc-point mesh. A plane-wave cutoff 
of 400eV is converged to 0.3meV/atom and the &-point mesh with Methfessel-Paxton smearing 
of 0.2eV is converged to lmeV/atom[ 10|. Projector augmented-wave (PAW) p seudopotential ifToTl 
calculations with the PBE generalized gradient approximation 17] give similar values, with a 
maximum error of O.leV (see supporting Table 1). From changes in supercell stresses for oxygen 
in different sites, we estimate the finite-size errors to be < 0.05eV; this is similar to the error 
found by using different computational cell sizes[10J. We use the climbing-image nudged elastic 
band lTT8l [T9l method with one intermediate image and constant cell shape to find the transition 
pathways and energy barriers between different interstitial sites. Along the path, the force is 
negated, while components perpendicular to the path are unchanged; the image relaxed to an 
extremum where the forces are less than 5meV/A, and restoring forces confirm that this extremum 
is a first-order saddle point. The attempt frequency prefactor for each transition is estimated with 
the Vineyard equation ll20ll . Only the restoring forces on the oxygen atom is used to compute the 
normal mode frequencies. We use the three vibrational modes from the initial state and the two 
highest modes at the transition state to estimate the prefactor. This does not include the coupling 
to and softening of neighboring Ti atoms; this underestimates the prefactors by less than 25% fl2TTl . 

Fig. [T] shows the hexagonal closed-packed unit cell of a-titanium and the three interstitial sites 
for oxygen. The crystal has space group 194, P6t. I mmc^22\. where the crystal basis d\ and a,2 are 
at an angle of 120° to each other in the hexagonal ("basal") plane with length = 2.933A, while 
the c axis is perpendicular to both with length c T ; = 4.63 8 A, and two titanium atoms per cell. The 
octahedral (o) site is the equilibrium configuration for oxygen and is surrounded by 6 titanium 
atoms in a symmetric arrangement, 2. 09 A away. The o- sites form a hexagonal lattice with a c- 
axis that is half of the titanium lattice. Atomic forces on oxygen in the unstable tetrahedral site 
displace it towards the basal plane, into the hexahedral site lfTOl ; the hexahedral (h) site is 5-fold 
coordinated and is 1 . 19eV higher in energy than the o-site. The three nearest titanium neighbors of 
the h-site are in the basal plane and are displaced to a distance of 1.92A, with two other neighbors 
directly above and below. The h-sites form another hexagonal closed-packed lattice as a-titanium 
with a translation of [000^]. The non-basal crowdion (c) site is 6-fold coordinated, but with 
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Site Wyckoff pos. [A] Z AE [eV] 

octahedral 2a (0,0,0) 2.09 6 +0.00 

hexahedral 2J(|,|,J) 192 5 +L19 

crowdion 6g(|,0,0) 2.00 6 +1.88 




FIG. 1. Wyckoff positions ll22l and relative energies for oxygen interstitial sites in a-titanium. The interstitial 
site energy is reported relative to the octahedral site energy. The geometry of each site is characterized by 
the nearest neighbor distance after relaxation, R nn , with the coordination number Z. Titanium atom sites are 
in white, while oxygen interstitial sites are in orange (octahedral), blue (hexahedral), and black (crowdion); 
the octahedral site is the ground state, with hexahedral and crowdion having site energies AE above. The 
octahedral sites make a simple hexagonal lattice; the hexahedral sites, a hexagonal closed-packed lattice; 
and the crowdion sites, a kagome lattice. 

lower symmetry and higher energy (1.88eV) than the o-site. The two titanium atoms that contain 
each c-site have been significantly displaced such that all six of the c-site's titanium neighbors 
are approximately 2.00A away. The c-sites form a kagome lattice Il23l in the basal plane and is 
repeated along the c-axis twice per titanium unit cell. The lowered symmetry of the c-sites in the 
kagome lattice mean that a distortion of the unit cell can give the different c-sites different energies. 
We also considered a crowdion site in the basal plane, which is unstable. A high formation energy 
is required to displace the two titanium atoms into the close-packed directions in the basal plane, 
while the two titanium neighbors of the non-basal crowdion can move in the softer pyramidal 
plane. Note also that oxygen retains its divalency in all three configurations, and in the transition 
states (see supporting Fig. 1). 
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FIG. 2. Oxygen interstitial sites and oxygen diffusion pathways in a-titanium. White spheres are titanium 
atoms, orange spheres are octahedral interstitial sites, smaller blue spheres are hexahedral sites, and the 
smallest black spheres are crowdion sites. The full transition pathway network is in the upper left, and is a 
superposition of all the remaining subnetworks. The individual transition networks are o<->o, o<-»h, o<->c, 
and h«->c. The bottom three networks are formed from pairing up transition networks o«->h, o«->c, and h<->c. 
All six networks contribute to the diffusion of oxygen. 

Fig.[2]show the interpenetrating network of transition pathways for oxygen between interstitial 
sites. There are two out-of-plane transitions (with rate /l 00 ) from each o-site with its two direct 
neighbors in the c-axis. The o-site is also surrounded by six h-sites and six c-sites and can transi- 
tion into them (with rates A h and A oc ). The h-site is surrounded by six o-sites and six c-sites and 
can transition into them (with rates Ah and A hc ). Each c-site resides in the center of a shared edge 
between two o-sites and two h-sites and can transition into them (with rates A co and A ch ). With the 
exception of the o<->o c-axis transition, all other pathways are heterogeneous, starting and ending 
at different site types, and have not been considered previously. The transition displacements and 
surrounding site symmetries are listing in Table [TJ 

Table [I] summarizes the symmetries and energetics of all possible transitions for oxygen dif- 
fusion in a-titanium. The transition rate from site i to j at temperature T is Arrhenius: Ajj = 
Vij exp(-Ejj/k B T), where E t j is the energy barrier and v, 7 is the attempt prefactor for the transition. 
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TABLE I. Transition pathways, prefactors v, and energy barriers E for oxygen diffusion in a-titanium, 
between octahedral (o), hexahedral (h), and crowdion (c) sites. The direction indicates the possible dis- 
placement vectors for the transition; the remaining transition vectors can be found by applying the point 
group symmetry operations. The symmetries are listed in Hermann-Mauguin notation: m is a mirror opera- 
tion through the basal plane (0001), I is inversion, 3 is a 3-fold rotation axis around the c-axis [0001] with 
inversion, 6 is a 6-fold rotation axis around [0001], and ^ is a 3-fold rotation axis around [0001] with mirror 
through (0001). The absolute rate of transitions h— >x and c— >x is thermally activated with the transition en- 
ergy barrier plus the site energy for h (+1.19eV) or c (+1.88eV), respectively; hence, all six heterogeneous 
transitions occur with similar absolute rates. 



Direction Symmetry v [THz] E [eV] 



o- 


->o 


(000 i) 


m 


11.76 


3.25 


o- 


^h 


<55°4-> 


3 


10.33 


2.04 


o- 


->c 


( uI 0) 

\663 U ' 


6 


16.84 


2.16 


h- 


->o 




3 
m 


5.58 


0.85 


h- 


->c 


(±0±±) 


2 

m 


10.27 


0.94 


c- 


^0 


( uI 0) 

\663 U ' 


1 


12.21 


0.28 


c- 


4h 


\6 U 64' 


I 


13.81 


0.24 



The barrier of the direct c-axis transition between o-sites — E QO — is too high to occur at relevant 
temperatures. However, the lower barrier o<->h transition also passes through a triangular face 
of three titanium atoms like the o<->o c-axis transition. This is similar to the instability of basal 
crowdion sites: the triangular face for the o<->o c-axis transition requires more energy to displace 
titanium atoms in the close-packed basal plane. The triangular face for the o<->h transition is in 
the softer pyramidal plane allowing for easier titanium atom displacement. Excluding the o<->o 
c-axis transition, all remaining transitions occur at approximately the same frequency — there is 
no single rate-controlling diffusion mechanism. As the probability of a site i being occupied is 
proportional to exp(-AEj/kBT) for site energy AEi, the absolute rate of transitions is proportional 
to exp(-(A£', + Eij)/k s T). The transition barriers from h- and c-sites are lower than from o-sites, 
but the occupancy probability for h- and c-sites are lower. Adding the site energy for h (+1.19eV) 
and c (+1.88eV) to the corresponding transition barriers reveal that all transitions occur with a 
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temperature dependence of about ~2.1eV; hence, all of the interpenetrating transition networks 
contribute to the diffusion of oxygen. 

We derive the exact rate equations for c-axis and basal diffusion using the multistate diffusion 
formalism [[24l l25l . The case of a single oxygen atom diffusing in a perfect a-titanium lattice is 
represented by the periodic unit cell in Fig. [T] with ten internal states. The connected network leads 
to basal and c-axis diffusion rates that are similar. The full diffusion equations (see supporting 
Eqn. Sll and SI 2) are simplified when the rates of escape from the hexahedral and crowdion sites 
are much faster than from the octahedral site; then, 



^basal — O, jj 



D r = c 



Ti 



3 1 /l h 

Ah + T^oc + -7- — Ac + O/loo 

4 4 Ao 

3 3 Ah 1 

-/loh + 0A oc + - —A. he + 7A0 
o Ao 4 



(1) 



(2) 



This assumes that the crowdion sites are able to fhermalize so that there is no correlated hops from 
the crowdion sites to neighboring sites. The contribution of the individual rates to diffusion is 
similar for A h, Ac, and Ac terms; at 300°C, the contributions are in ratios of 13.3: 1.45: 1 for basal 
diffusion, and 3.3:0:1 for c-axis; at 600°C, 7.1:1.8:1, and 1.8:0:1; at 900°C, 5.3:2.0:1, and 1.3:0:1; 
and at 1200°C, 4.4:2.1:1, and 1.1:0:1. In all cases, the rate A Q0 is significantly smaller, contributing 
only < 10~ 4 at 1200°C. Over the temperature range of interest, all of the heterogeneous networks 
contribute to the diffusion of oxygen. 

As temperature increases, we do not expect the crowdion site to achieve thermal equilibrium; 
hence, an increasing fraction of the o— >c jumps will become correlated basal o— >o jumps, and h— >c 
jumps will become correlated h— >h jumps. This high temperature behavior can be approximated 
by removing the crowdions as metastable states from the network, and using A oc as the rate for 
direct basal o— >o transitions (and similarly /l hc for h— >h transitions). Then, the diffusion rates are 
bounded above by 

. r 1 1 a., 1 

(3) 



3 1 /i Q h 

Ah + 7 Ac + 77 — ^hc + 0Ao 
I I /(ho 



D highT 2 
« c Ti 



3 3 /i h 1 

-/loh + 0/loc + T— Ac + 7A0 
O 4 Afro 4 



(4) 



At 1200°C, the high temperature Eqn. [3] is 41% larger than Eqn. [T] and Eqn. [4] is 48% larger than 
Eqn.[2| at 300°C, the differences are only 16% and 23%. This suggests a small underestimation of 
diffusion rates at the highest temperatures. 

Fig. [3] shows the diffusion coefficient from the multistate diffusion equations against experi- 
mental oxygen diffusion data. From the diffusion equations Eqn. [T] and Eqn. [2j the temperature 




1000 / T [K" 1 ] 

FIG. 3. Analytical results and experimental data of oxygen diffusivity in a-titanium. We compare our 
analytical DFT model (bold line) to experimental data from the literature survey by Liu and Welsch[26] 
(thin lines), and experiments by Bregolin[27 ] and Vykhodets l|28|| (symbols). Over the temperature range of 
300-1200°C, the diffusion rate is Arrhenius D = 2.18xl(T 6 m 2 s -1 with E act = 2.08eV. 

behavior should follow the barriers, E h, E oc , and E hc + E ^ - E^o- An Arrhenius model with D Q 
= 2.18 x 10~ 6 m 2 /s and a single barrier E act = 2.08eV matches Eqn. [I] and Eqn. |2] to within 15% 
over the range 300-1200°C, with largest deviations at low temperatures. The experimental data 
come from the literature survey lEoll by Liu and Welsch and more recent experiments Il27l 1281 using 
nuclear reaction analysis. The activation energy matches well to experiment while the absolute dif- 
fusion coefficient is a factor of ten below experimental values — well within the expected accuracy 
of density-functional theory for diffusion. 

First-principles calculations determine the pathways for oxygen diffusion in or-titanium, includ- 
ing a new interstitial non-basal crowdion site for oxygen in titanium. Other than the high-barrier 
direct c-axis transition between octahedral sites, all transition paths are heterogeneous (o<->h, o<->c, 
and h<->c) and contribute to diffusion over a wide temperature range. This shows that even well- 
studied materials science problems can have surprises: new configurations and new transitions 
give rise to complexity for single atom diffusion. Moreover, the new sites suggest interesting in- 
teractions with titanium vacancies are possible, as they should destabilize nearby crowdion and 
perhaps hexahedral sites. We expect other interstitial elements like carbon and nitrogen to have 
similar diffusion networks in titanium, and in other hexagonal-closed packed metals like magne- 
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sium and zirconium. This new understanding of oxygen in titanium can serve as the basis for 
controlling oxygen diffusion in alloys, growth of oxide phases in titanium, and related challenges. 

This research was supported by NSF/CMMI CAREER award 0846624 and Boeing. The au- 
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